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Abstract 



We report on a high temperature perturbation expansion study of the 

superfluid-density spatial correlation function of a Ginzburg-Landau-model 

superconducting film in a magnetic field. We have derived a closed form which 

expresses the contribution to the correlation function from each graph of the 

perturbation theory in terms of the number of Euler paths around appropriate 

subgraphs. We have enumerated all graphs appearing out to 10-th order in the 

expansion and have evaluated their contributions to the correlation function. 

Low temperature correlation functions, obtained using Pade approximants, 

are in good agreement with Monte Carlo simulation results and show that 

the vortex-liquid becomes strongly correlated at temperatures well above the 

vortex solidification temperature. 
PACS numbers: 74.40.-hk, 74.60.-w 
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I. INTRODUCTION 



Because of the combination of high transition temperatures, strong anisotropy, and short 
coherence lengths which occurs in high temperature superconductors, strong thermal fluctu- 
ations are present over a wide temperature interval in these materials. Thermal fluctuations 
are especially important in a magnetic field where they are responsible for the melting of the 
AbrikosovB vortex lattice at temperatures below the mean-field critical temperature giving 
rise to a vortex liquid stateS. In this paper we report on a study of correlations in the vortex 
liquid in the extreme anisotropy limit of decoupled layers. We have evaluated the lead- 
ing terms in the high temperature perturbation expansion of the superfluid density spatial 
correlation function for the Ginzburg-Landau model of a superconductor film. 

The thermodynamics of this system is unusual because of Landau quantization of the 
order parameter fluctuations!; our calculations are carried out within the lowest Landau 
level approximation in which only the lowest gradient energy fluctuations are retained. This 
approximation is valid near the mean-field transition temperature and is ordinarily valid 
throughout the vortex liquid state, although fluctuation effects may be strong enough in some 
high temperature superconductors to drive the solidification transition outside of its range of 
validity. High-temperature perturbation expansions for the free energy of this model, even 
when evaluated to high ordei0i, exhibit little evidence of the Abrikosov vortex lattice state 
which is expected to occur at low temperatures. Recent Monte-Carlo simulationsiM, on the 
other hand, generally obtain results indicative of a weak first order phase transition between 
the two states. (See however Ref. [ |^ where a different conclusion is reached.) We find 
that the perturbation expansion for the superfluid density correlation function, unlike that 
for the average superfluid density studied by previous workers!'!, provides clear evidence of 
strong correlations in the vortex liquid which presage the appearance of an ordered state at 
low temperatures. 

The paper is organized as following. In Section II we introduce the quantity we study, 
(|A(g)p) which is proportional to the Fourier transform of the superfluid-density spatial 
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correlation function. In Section III we outline the high temperature perturbation expansion 
for this quantity and present our closed form result for the contribution to it from individual 
diagrams which appear in the expansion. In Section IV we discuss our evaluation of all terms 
out to tenth order in this expansion and discuss some comparisons with expansions for the 
free energy performed by earlier workers which allow us to check our results. Extrapolations 
of our finite order results to low temperatures using Pade approximants are presented and 
compared with Monte Carlo simulations. Section V contains a brief summary. 



II. SUPERFLUID DENSITY CORRELATION FUNCTION 

The free energy per unit area of the Ginzburg-Landau model superconducting film in a 
perpendicular magnetic field is given by 

/[vl>] = airm' + f + ^l(-^^V - 2eAm' (1) 

where \l/(r) is the order parameter. \l/(r) is proportional to the wavefunction for the cen- 
ter of mass of the Cooper pairs0 and 5 = V x A (We employ the Landau gauge so 
that A = {0, Bx,0).) The quadratic terms in Eq. (|ip are minimized by order-parameters 
which correspond to a lowest Landau level (LLL) wavefunction for the Cooper pairs. In 
using Eq. (|l]) we are neglecting screening so that B is the external magnetic field. This 
approximation is valid for type-II superconductors near the upper critical field, the regime 
of interest in this article. The mean-field-theory superconducting instability occurs when the 
quadratic terms in Eq. (||) become negative for the LLL, i.e. at T^^^ where aniT^^^) = 0. 
[an = a + heB/m* . ) In the LLL approximation we assume that fluctuations in higher Lan- 
dau level channels can be neglected or at least absorbed in a renormalization of a 
Then the order parameter \l/(r) can be expanded as 

m = E C,(L,)-V2(^)i/4eWe-¥(-^)^ (2) 

where the number of terms in the sum over p is A^,^ = LxLy/2iTi'^ = L^Ly^eB / -kK) . This 
expansion leads to the following expression for the Ginzburg-Landau model free energy: 



P P1P2P3P4 i=l i=l 

XCpiC'p2^P3^P4'^Pl+P2,P3+P4} (3) 

where an = C({T){1 — H/Hc2(T)). Fluctuation effects in the model are regulated by the 
dimensionless parameter, g = aniT^^"^ / PksTY^'^ . 

The central quantity in our work is the superffuid-density spatial correlation function 
which we define by 

Xsfd{R) = {\^{m^{r + R)\') - {\^{f^\'){\^ir + R)\') (4) 

Translational invariance of the system guarantees that the right hand side of Eq. (|) is 
independent of r. The modulation of the average superfluid density when the homogeneous 
system is disturbed by weak pinning can be expressed in terms of xsfd{R)- 

S{\^{^\') = ^J d^XsFDi^ - f)6a{r') (5) 

where 6a{r') reflects the modulation of the mean-field transition temperature by pinning. 
The Fourier transform of XsFD{r) is 

XsFDiq) f d'rf d'r'mif^l'l^if'l') - {\^ {f\') {\^ {f')\')] exp[^q ■ (r-r^)] (6) 

LxLy J J 

We evaluate Xsfd{(1) by expressing it in terms of 

= ^Y1 Cp,Cp^Sp^^p,+g^ exp[-ifq.x{pi + p2)/2] (7) 

SO that 

XsFDiq) = ^[exp[-g^£V2](|A(g-)H - (A(g = 0))^] (8) 

A(g) satisfies the following informative sum rulei for each configuration of the Ginzburg- 
Landau system, 

^E[iA(0r-i/iV0]=o (9) 
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where A{q) = A(g)/Ao and Aq = A(g = 0) is proportional to the integrated superfluid 
densitylii. (Note that A(g) is invariant when the order parameter is multiplied by an overall 
constant.) (|A(g)p) is a particularly revealing quantity to examine in studying correlations 
in the vortex liquid. Eq. (|) guarantees that (for large N^) limg^oo N(i,\A{q)\'^ = 1 for 
any vortex liquid configuration. For example, in the high temperature (vortex gas) limit 
(|A(g)p) = for all q ^ 0. On the other hand it is easy to show that in the low 

temperature limit the mean-field Abrikosov lattice configuration of the order parameter 
gives 



\m\' 



1 aX q = G, 
otherwise. 



(10) 



where G is any reciprocal lattice vector. Thus N^{\A{q)\'^) = sv{q) shows exactly the 
behavior which would be expected for the static structure function of a classical fluid with 
particles in both low temperature (solid) and high temperature (gas) limits. It seems 
clear that sy(g) must be closely related to the static structure factor of the zeroes of the 
order parameter, the vortices, although we do not believe that they are identical at all 
temperatures. For the purposes of the present study it is sufficient to observe that hv{q) = 
sv{q) — 1, which we call the vortex correlation function, is a convenient measure of of the 
degree of correlation in this system. 

III. HIGH TEMPERATURE PERTURBATION EXPANSION 

|A(g)p can be expressed in the form 

\Aiq)\' = N^' Y: Cp,C,,C,3C,,5,,+,„,3+,,(5,,,,+,^e'^''^^(P-^^)). (11) 

P1P2P3P4 

At high temperatures we can evaluate its thermal average 

X{dC,dC,\m?e-^^''^^^' X e-^/*^l*l' (12) 
where Z is the partition function 



Z= j \{dC,dC,e-^ ^ '"^^^^ X e-^^"^^^*!' (13) 

by expanding the contribution to thermal weighting factors from the quartic contribution 
to the Landau-Ginzburg free energy. The perturbation series can most easily be handled in 
terms of Feynman diagram At n-th order , there are contributions from diagrams with 
(n + 1) vertices and 2{n + 1) edges in which the edges represent the Gaussian approximation 
correlation functions, n vertices represent |^|^ terms proportional to (3 and the additional 
'external' vertex corresponds to |A(g)p. An important consequence of the Landau quantiza- 
tion of order parameter fluctuations is the fact that the Gaussian approximation correlation 
functions, 

{C,,C,) = 6p,,pkBT/aH, (14) 

are independent of the momentum p. Each vertex has two directed outgoing lines to represent 
the factors of Cp associated with it and two incoming lines to represent the factors of Cp. All 
diagrams containing single-loop dressings of edges or vertices (except the external vertex) 
can be eliminated by following Ruggeri and ThoulesJ^ and expanding in terms of 'Hartree- 
Fock' approximation correlation functions. This renormalization replaces^ an in Eq. (0) by 
a = —2aH/{g'^±[g^+'ig'^Y^'^). (The + (— ) sign applies for (7 < {g > 0).) a remains positive 
for all temperatures so that the renormalized expansion parameter of the high-temperature 
perturbation expansion, x = {fiksT) / [AT^l'^a^), remains finite at all temperatures, (a ~ uh 
for 5( > 0, a ^ \aH\/g'^ for ^ 0; a = a//(l — 4x)"\ g"^ = {Ax — l)^/4x and x = 1/4 at 
jiMF^ When the expansion is performed in terms of the Hartree-Fock correlation functions 
(|A(g)p) is the sum of all the connected diagrams without single-loop dressings except at 
the external vertex. 

It is convenient to label the external vertex as vertex 1 and to label the incoming momenta 
at vertex i as p2i and P2i-i- Then the contribution to (|A(g)p) from an n-th order diagram 
IS given by (fcBT/a)2(-x)"A/'^i/(g)/n! where 

I{q} ^ e-^'''?^(^^+^'?^)(f )-"/2 Ep, / dps... J dp2n+2 W^tl S{J:^{M,, - N^,)p,) 



X (V.-f.. + ^p^+.Je-^^fe+^''^)^^ exp{-f E';±liP2, - P2,-i?}, (15) 

p[ and P2 are the two outgoing momentum labels at vertex 1, /i labels the vertices, and 
M^i is unity if i goes into the vertex fi and is zero otherwise, while N^i is unity if i comes 
out of /X and is zero otherwise.! To obtain this result we have noted that the integral is 
invariant under a shift of all momenta, set pi = and multiplied by A^^. It turns out that 
the integral I{q) can be evaluated exactly and expressed in terms of the number of Euler 
paths in the two subgraphs obtained by deleting the external vertex and making the possible 
contractions. The two contractions correspond to the two delta functions which fix either 
p'l or P2 in Eq. (|T5|). We denote the Euler path numbers for the two subgraphs by and 
T^. The following result is derived in the Appendix: 

m = {^^M-^<l'n + ^exp(-^g2£2)} (16) 
if T^xT^ ^ 0, and 

/(g-) = (iVAlo + l)/(T^ + T^) (17) 

if T^xT^ = 0. 

Using the above result for I{q) and writing the number of appearances of a given graph 
as 4"+^?7,!/G„_|_i g we obtain the following formally exact expression for (|A(g)p): 

= — (^)2(1 + N^6,^,o - 4a:exp[-^] + exp[-^] 

+8x='(l + Ar<^%o) + 16x^(0.5 exp[-^] + expf-g^^^]) ^ (^g) 

where g labels graphs and the sum at n-th order is over (n + l)-vertex graphs. an,g{q) is 
given by 

= \ " (19) 

Gn+i,,(i;<i,+r„^,)(^ + A^^^g-o) If T4 x T^^ = 0, 



where T^^ and T^^ are the number of the Euler paths of the two contracted n- vertex graphs 
and Gn+i,g is the number of automorphisms of the {n + l)-vertex graph with one external 
vertex. The diagrams which appear up to second order in the series and their associated 
properties are hsted in Table |. The explicit expression in Eq. (|T8]) can be confirmed from 
the entries in this Table. 

We observe in Eq. ([19|) that contributions which survive to the large |g| limit come only 
from graphs where T^^ equals zero or T^g equals zero and that the only terms in (|A(g)p) 
which are independent of A'"^ come from the same set of diagrams. The contribution inde- 
pendent of N(f, is The remaining contributions to (|A(g)p), which are proportional 
to A*"^^, contribute to Xsfd{.<1) and are due to correlations in the thermally fluctuating 
superfiuid density. We thus obtain explicitly from the perturbation expansion that 

lim(|A(g-)r) = (A^)/iV,. (20) 

This result was claimed earlier on the basis of the sum rule (Eq. (^). 
The real space correlation function is given by 

(|vl/(f)nx[/(f+/?)p) = ^/ ^/g(|A(g-)r)exp[-g2f/2]expHg-/?] (21) 

= {^n^fm:<^nAR){-^^T: (22) 

where 

anAR) = r T (^M-^^ijY] + exp[--S^(f )^]). (23) 

In Eq.(p3D Tn+i,g = T^^g + T^g is the number of Euler paths in the uncontracted {n + 1)- 
vertex graph. We see again here that only graphs with T^g = or T^g = remain finite for 
R ^ oo where correlations vanish. The contribution of this subset of graphs is (|\E'(r)p)^. 
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IV. FINITE ORDER RESULTS AND EXTRAPOLATION TO LOW 

TEMPERATURES 



We have written a computer program which generates all relevant graphs represented 
by their adjacency matrices.0 From the adjacency matrices we calculate the number of 
Euler paths for the subgraphs, (T^^ and T^c,) and from the graph generating algorithm 
we calculate the symmetry factor Gn+i,g- In Fig. (|lD, we show one of the graphs which 
appears at second order in the expansion, the two graphs which result from the deletion 
of its external vertex, and the adjacency matrices of all three graphs. (The number of 
Euler paths equals the determinant of the minor of matrix appropriately formed from the 
corresponding adjacency matrix. 0) In this way, we have evaluated the series exactly up to 
tenth ordei0. We have checked our results for (|A(g)p) by confirming that the sum rule 
Eq. (H), is satisfied and that the results for both (|Aop) and I]g |A(g)p exp{ — ^g^f^} are 
correct order by order in perturbation theory. The latter two quantities can be related to 
derivatives of the free energy. The free energy is given by 

F = -ksTln Z = kBTN^{\n{a/7TkBT) + hn^x)) (24) 

where the perturbation expansion for f2D{,x) was first calculated by Ruggeri and Thoules^l 
and the expansion was extended to order x^^ by Brezin, Fujita and Hikami.i Differentiating 
Eq. ( P^D once with respect to [3 and find that 

E(|A(g-)r)exp{-ig^/n = {^f x ^ + (25) 

Similarly differentiating twice with respect to an gives 

/|A |2\ /A \2 1 (kBT ^ {I - Ax){l - 2xf^o (x) ) 4x(/^g(x) 

^ " ^^'^ - W/^^ ^ OTi^? ^ ^ ^ 

where (Aq) = (A;bT/5)(1 - 2a;/^^(a:))/(l + Ax). 

The asymptotic high-temperature expansion of (|A(g)p) can be extrapolated to low tem- 
peratures using Fade approximants to describe the x dependence at each value of q. Com- 
parisons between Fade approximants of the series and our Monte Carlo simulation resultJl 



are given in Fig. ^for different temperatures. shows a well-defined peak for q near |G|, 

where G is a reciprocal lattice vectors of the vortex solid, Such a peak is characteristic of a 
strongly correlated liquid. At low temperatures the [2n, 2n] Fade approximants appear to 
overestimate the correlations while the [2n+l,2n+l] approximants appears to underestimate 
the correlations. Quantitative agreement with Monte Carlo simulation results is obtained 
for temperatures above T^^^ and for temperatures below T^^^ with g"^ < 'b. A continuous 
phase-transition to a vortex solid state at low temperatures would be indicated by the di- 
vergence of Sy{q) at q near |G| and hence by a pole in the Fade approximant to (|A(g)p) 
along the positive real x axis. It might be significant that for the [4,4] Fade approximant 
poles on the positive real axis appear over a finite range of q near \G\. In Fig. ^ we plot 
the q dependence of the value of Xc at which these poles occur. The highest temperature 
at which a pole occurs is x = 9.69 corresponding to g"^ = 36.8 for q = 0.91|G'|. This value 
of g^ compares with the value g^ = 43.5 ± 1.0 at which a weak first order phase transition 
occurs according to the Monte Carlo calculations.! No poles on the positive real axis for x 
occur in the [5, 5] Fade but we expect that they will recur at slightly lower temperatures in 
the [6,6] Fade. We believe that the structure in sy(g) is a precursor of the solidification of 
the vortex lattice at low temperatures. The proximity of the vortex lattice state is apparent 
in the perturbation theory for Sy(g) but is hidden in the perturbation theory for the free 
energy because of the weakness of the thermodynamic singularity associated with the phase 
transition. However, accurate estimates of the transition temperature based on the pertur- 
bation expansion for sy(g) would require calculations to be carried out to higher order than 
we have found possible to date. 



V. SUMMARY AND CONCLUSIONS 

We have studied the superfluid density spatial correlation function (|A(g)p) using high 
temperature perturbation expansions and extrapolated our results to temperatures below 
the mean-field transition temperature by means of Fade approximants. Good agreements 
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with Monte Carlo simulation data is obtained for g > — -^/S. Our results demonstrate that 
the vortex liquid is strongly correlated below T*^^. A result argued for previously on the 
basis of a sum rule which related the large wavevector limit of the correlation function 
to the average superfluid density was obtained explicitly from the perturbation expansion. 
We argue that these perturbation expansion studies bespeak the phase transition to a two- 
dimensional vortex solid that is believed to occur at lower temperatures. However, higher 
order calculations than we have been able to complete to date would be necessary in order 
to obtain estimates of the transition temperature which are competitive in accuracy and 
reliability with those obtained previously from Monte Carlo simulations of the same model. 
This work was supported by the Midwest Superconductivity Consortium through D.O.E. 
grant no. DE-FG-02-90ER45427. The authors are grateful to Steve Girvin, Brendan McKay 
and Lian Zheng for helpful interactions. 



APPENDIX: 

The integral to be evaluated has the form: 

I{q) ^ e-'^^^(^^+^^^)(|f Ep, / dp,... J dp2n+2 n;:l2 S{J:^iM,, - N,,)p,) 

X + S,'^+Je''''^^^^''^^P^ exp{-f E;tl{P2, - P2,-ir}, (Al) 

Of the (2n + l) variables pi only n are independent because of the {n + 1) delta functions. 
It is convenientlll to choose = P2fj.—p2fi-i (yU = 2, ■ ■ ■ , n + 1) as the n independent variables. 
We label the contributions to I{q) from the two choices for the delta function at the external 
vertex as I^{q) and I^{q) respectively. We use this delta function to eliminate the sum over 
P2 = Si which can be expressed as a linear function of the independent variables. For X =A 
or B 

n+1 
fj.=2 

The integral over the independent variables can then be expressed in terms of these coeffi- 
cients. The Jacobian for changing variables from p^ to s^>i is independent of qy] at = 
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the change of variables is identical to that required for the diagram obtained by deleting the 
external vertex and contracting the two outgoing edges at that vertex with the two incoming 
edges. (The two ways of doing the contraction correspond to I^{q) and I^{q) respectively.) 
Following the work of McCauley and ThoulesJli we note that the inverses of the required 
Jacobians equal the number of Euler pathli for the n-vertex graphs which results from the 
deletion of the external vertex and the two possible contractions. We denote the Euler path 
numbers by . The integral over the is then elementary and we obtain 



Where nf = -1 + - E;±l{t^y, = -I + - \ E^Uit^)' and = i E^Uit^f. 
However, we know that I^{q) is real which leads to the requirement that = and implies 
that I^iq) = exp(— — l)/2)/T^. The value of can be inferred by noting that 
exp(— g^^^/2)/^(g) is proportional to an integral which appears in the expansion of the 
free energy and equals N^/T where T is the number of Euler paths in the original graph 
before deletion. It follows that a-^ = T/T^ and hence that 

H<r) = {^^M-^q'f") + ^exp{-^q'f)}. (A4) 

(Note that T^+T^ = T.) This analysis fails in the special case where x = 0. For that 
case we have I^{q) = n^^^fi or I^{q) = and a similar analysis gives the corresponding 
= N^/T or = 1/T and so 

/(g) = (iV^%,o + l)/(T^ + T^). (A5) 
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TABLES 

TABLE I. All diagrams up to second order in the high-temperature perturbation expansion. 

The open circle in each diagram represents the external vertex. 
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FIGURES 

FIG. 1. Example of the adjacency matrix representation of a graph which appears at second 
order in the expansion and of tlie ealeulation of Euler path numbers of the grapli. The graph is 
represented by an adjacency matrix A with elements a^- where aij is equal to the number of edges 
going from vertex i to vertex j. The number of Euler paths is equal to the determinant of the 
minor of the matrix B with elements defined by bij = Sij^i^Uik — aij. The open circle in the 
uncontracted graphs represents the external vertex. The two graphs on the left result from the two 
possible contractions after deletion of the external vertex. The contribution of this graph to the 
series is 02(5) = 0.5exp(-g^^^/4) + exp(-g^£^). 

FIG. 2. Comparisons of the perturbation calculation and Monte Carlo simulations. The y-axis 
label is (|A(g)p), while the x-axis is the wave vector in the units of the inverse magnetic length 
The upper- left panel is for a temperature above Tj^^, while the other panels are for temperatures 
below T^^ . The arrows in the plots indicate the location of reciprocal vector, |(5|, of the vortex 
lattice state. The Fade approximants sometimes behave poorly at small wavevectors. {q(, ^ 1.0.) 
We have 'clipped' anomalous small wavevector behavior in some of these figures so as not to obscure 
the behavior for q near G. 

FIG. 3. Wavevector dependence of one of the poles in the [4, 4] Fade approximant. The pole 
moves from negative to positive x by passing through 00. This approximant has poles for q near 
G at temperatures which are somewhat higher than the temperature where a (weakly) first order 
phase transition is indicated by Monte Carlo simulations. (See text.) The arrows in the plot 
indicate the location of reciprocal vector |G| of the vortex lattice state. 
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